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Abstract. We show how to descritize the Korteweg-de Vries (KdV) equation 
in such a way as to preserve all the Lie point symmetries of the continuous 
differential equation. It is shown that, for a centered implicit scheme, there are 
at least two possible ways of doing so. 



1. Introduction 

Symmetries play an important role in our understanding of nature, they are intrinsic and fun- 
damental features of differential equations in mathematical physics. Hence, they should be 
retained when discrete analogs of such equations are constructed. The aim of this paper is to 
give a systematic approach for discretizing partial differential equations (PDEs), involving one 
function of two variables u(x,t), such that all the Lie point symmetries of the original problem 
are preserved. To illustrate the procedure we investigate the case of the KdV equation. We 
will show how it is possible to obtain invariant implicit schemes. Explicit schemes can also be 
derived but they often exhibit numerical instability. Finally, a numerical application is realized 
using the schemes obtained. 



2. Discretization procedure 

The discretization of a PDE involving one depend variable u and two independant ones, x 
and t, consists of sampling, couples (x,t) in the real plane. Each point is labelled by a set of 
discrete indices: 

(Xm,ni tm,n 

), m,nEZ. (1) 

The discretization of the independant variables induces a natural discretization of the depen- 
dant one 

Um,n M(Xm,mtm,ii) ■ (2) 
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A two-variable partial difference scheme (PAS) approximating the original PDE will be 
a set of three equations relating the quantities {x,t, u] at a finite number of points 

Ek\\Xm+i,n+jjtm+i,n+jiU m +i jn -\-jy) = 0, 

\<k<3, h<i<i2, j\<j<h, 

such that, in the continous limit, all three equations © reduce to the original PDE. By choice, 
we suppose that E12 = describe the lattice and impose that in the continous limit these 
equations go to the identity = while £3 = goes to the PDE. We also require that the 
system © be invariant under the same group of symmetries as the original PDE. 



2.1. Lie point symmetries of PAS 

Let G be the group of Lie point symmetries, of order N, for a given PDE 

=Xi(x,t,u,e) =x + eE, l (x,t,u)+0(£ 2 ) 
t\ = Ti(x,t,u,e) = t + ex 1 (x,t ,u) + 0(e 2 ) , 
u] = Ui(x,t,u,e) = u + e^ l (x,t,u)+0(e 2 ). 

To this group is associated the Lie algebra of vector fields 

V/ = % l (x,t,u)d x + X 1 (x,t,u)d t + § l (x,t ,u)d u , 



1 = 1,. ..,N 



(4) 



1 = 1,.. .,N. 



(5) 



As in the continous case, we define a prolongation of ©, but in a different fashion. The 
prolongation is realized by requiring that the vector field acts at all points figuring in ©: 



where 



m+12 n+jl 

P rV,:= £ £ ^ J d Xij + xj J d tij + ^d Uij \ 



^ij = ^\ x ij^ij^ u ij)^ r/ y r'f.v, ,.//,.//,/) and $1 = (j> l (xij,tij,Uij) 



(6) 



Let M be the manifold on which G acts, i.e. M ~ {x m+i] , t m+h )W+ y, , u m+h )W+ y, , . . . , 
Xm+i2,n+ ]2 •> ^m+;'2,n+j2 ' ^ m +i 2 ,ii+j 2 }- The quantity / : M — > E is said to be strongly invariant if it 
satisfies 

prV / / = 0, l = l,...,N. (7) 

Using the method of the characteristics, we obtain a set of elementary invariants {I\,...,I a }. 
The number of them is given by 



a = dim M — rank Z, 



a>0, 



(8) 



where Z is the N x (z 2 — h + h~ h ) matrix 
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+ I 2, n +J2 m — ^ ; 2!»+72 



d) 1 ■ \ 

rm+i2,n+ /2 



§m+i2-n+h) 



(9) 



?m+i2,n+j2 m+n,n+]2 



A difference equation will be strongly invariant under G if it can be written as 

E{h,...,I a ) = 0. 
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(10) 



Other invariant difference equations can be obtain if there exist expressions 

E {Xm+ii ,n+ j\ ; tm-\-i\ ,n+ 71 5 ^m+ii ,n+ j\i • • • )-^m+/2,«+j2 ' ^wi+J2,'!+ j'2 ' ^w+!2,n+ 72) = 

such that the rank of Z is not maximal. Such equations are said to be weakly invariant and 
satisfy 



prV/£ 



0, / = !,...,#. 



(11) 



E=0 



3. Discretization of the Korteweg-de Vries equation 

It is well known [2|, that the KdV equation 

u t = uu x + u xxx (12) 

admits a 4-parameter Lie symmetry group of point transformations, generated by the infinites- 
imal operators: 

Vi = a*, V 2 = 3, V3 = td x -d u , 
V4 = x<9 x + 3td t — 2ud u . 

Before computing the strong invariants of ( TOl . we point out the fact that 

7+ = tm,n+l ttn,n (14) 

is a weakly invariant equation. Hence, the strong invariants can be computed on a scheme with 
flat time layers, which is physically desirable. For numerical stability reasons, we compute the 
invariants on an implicit scheme involving the points shown below 



(13) 



(x ,t,u ) (x~,t, U-) 



(x,t,u) (x + ,t,u + ) (x ++ ,t,u ++ ) 
• • 




For notational simplicity, we write 

(Xm,nitm,niMm,ji) — [X,t,U), (<Xm+l,n>^m+l,n> ^m+l,n) = \^i^^U)^ [Xm,n±l 5^m,)i±lj^m,«±l) = C*±> ^±) 

and so on. The results of the computations are: 

h + h ++ h- h + 
h = j— n = -j — ^3 = j — h = — 
h- h + h t + 



T o + tu T tit* \ 

h = —? I 6 = h + {u + -u) 

h+ 

I U— —U \ I U—U. 



n+ ) \ h. 
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(15) 



From the set of invariants (fT31) . different type of schemes can be obtained. We will 
restrict our attention to schemes for which the evolution of the solution is given implicitly over 
a lattice that can be determined explicitly. 

3.1. Lattice with uniform steps in x 



T+ = 
h = \ 

h = \ {h + hh) {b +h) + (ho - b) - (h - h) ■ 

In terms of the original variables, system (fToT) reduces to 

T+ = 

h + = h- = h 

/a \ / a a \ / a rs a, . rs a, a \ / a a ^ 

(u — u) — U-) — 2w+ + 2u- —u ) a [u+ — U- 

— u ~ I at r 1 



(16) 



(17) 



t 2h 2h 3 t 2h 

The third equation of (fTTt converges in the continuous limit to the KdV equation if 

<T^0 



\R\ < oo. 



(18) 



(19) 



Hence <7 = O(t). Further more, the equations for the lattice can be solved and give 

t m ,n = Y(m), 
x m ,n = a(m)n + P(m), 

where y(m), a(m) and /3(m) are arbitrary functions that we can choose. The extra term 

O (w_|_ — U-) 

in the discrete KdV equation is geometrically seen to be a correction term to the solution taking 
in account the possible displacement of the points in x as the time evolves. 



3.2. Lattice depending on the solution 

Schemes for which the evolution of the lattice depends on the solution have already been 
obtained by Dorodnitsyn et al , HI. They use the concept of Lagrangian coordinates which is 
a point of view that differs from ours. Nevertheless, we get 



(20) 



h = 



ho ~h h~ h 



/ 2 + 2 + 2/f 1 + (/ 3 / 1 )- 1 l/ 2 + l 



h 
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Figure 1. Numerical solution on the evo- 
lutive lattice scheme. 



Figure 2. Evolution of the lattice. 



which in terms of the original variables gives 

T + 
o = 

u — u 6 

T h ++ + 2h + + 2h_+h 

From the numerical point of view, the scheme d2TTl is interesting since there is no nonlin- 
ear term in the equation governing the evolution of the solution. Hence, the computation time 
is reduced since we avoid nonlinear numerical methods. The absence of the nonlinear term is 
compensated by the specific evolution of the lattice. 





-xu 



(21) 



}++. 



u x -u x 



[h+++h-\. h + /z_ 



4. Numerical application 

This section is not meant to be an exhaustive investigation of all the numerical aspects of our 
invariant schemes. Actually, the aim is to give an illustration of the usefulness of such schemes 
for computing numerical solutions. 

The solution we have choosen to investigate is 

u(x,t) = — . (22) 

Now, if we want to use the scheme (fPTt . we must first specify the lattice over which the 
numerical solution is to be computed. From all the possible choices, we will look at two 
specific cases. As a first choice, since d22l is linear in x, we may take 



The solution of d23l is given by 



x 

X-m+l,n Xm,n = C = X— . (23) 



t + x\ m 

h I -—} n (24) 



where ho is the initial step in x. With this choice, scheme (fT7t reduces to the scheme d2"TT) . 
As a second choice, we choose o to be zero. In that case, the lattice is orthogonal and the 
discretization reduces to the standard implicit discretization of the KdV equation. 
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„ -5 
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UJ, 
o 

~~ -10 




Figure 3. Evolution of the error as a 

fonction of z: h(to) =0.1, evolutive 

lattice, orthogonal lattice. 
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Figure 4. Evolution of the error 
as a fonction of h : z = 0.000001 

evolutive lattice, orthogonal 

lattice. 



In the Figures [T] and |2] we have plotted the evolution of the solution and lattice for the 
scheme define on the evolutive grid. The initial step in x was ho — 0. 1 and in the time variable, 
we've taken a constant step, x = 0. 1 . 

We've also plotted the evolution of the error as a fonction of the step sizes in the figures 
|3] and IU for the two schemes choosen. The error has been computed using the /«, norm. As it 
can be seen, the error due to the evolutive lattice is much smaller than the one generated with 
the orthogonal mesh. 

From this numerical simulation, we clearly see, that for the scheme (fT7t . the choice of 
the lattice influences a lot the numerical precision. At the moment, we don't have a systematic 
way of determining the lattice that will give the best numerical results. On the other hand, with 
the scheme (OTT) everything is determined by the system. The only unknowns are the initial 
steps in x and those in t. 



5. Conclusion 



We have shown that it is possible to discretize the KdV on an implicit scheme such as to 
preserve all the symmetries of the original equation in at least two ways. From the theoretical 
point of view, these schemes are interesting since they preserve all the group properties of the 
original equation. Hence, symmetry reduction can hopefully be used to obtain exact solutions. 
From the numerical point of view, we see from our exemple, that the invariant schemes have 
the potential of giving better numerical results. 
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